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Abstract 

In this paper we use a similarity transformation connecting some families of Nonlin- 
ear Schrodinger equations with time- varying coefficients with the autonomous cubic 
nonlinear Schrodinger equation. This transformation allows one to apply all known 
results for that equation to the original non-autonomous case with the additional 
dynamics introduced by the transformation itself. In particular, using stationary 
solutions of the autonomous nonlinear Schrodinger equation we can construct ex- 
act breathing solutions to multidimensional non-autonomous nonlinear Schrodinger 
equations. An application is given in which we explicitly construct time dependent 
coefficients leading to solutions displaying weak collapse in three-dimensional scenar- 
ios. Our results can find physical applicability in mean field models of Bose-Einstein 
condensates and in the field of dispersion-managed optical systems. 



Key words: Nonlinear Schrodinger equations, similarity transformations 



Preprint submitted to Elsevier Science 



8 February 2008 



1 Introduction 



Nonlinear Schrodinger (NLS) equations appear in a great array of contexts 
[1] as for example in optics of nonlinear media [2], photonics [3], plasmas [4], 
mean-field theory of Bose-Einstein condensates [5] , condensed matter physics, 
[6] and many other fields. In some of these fields and many others, the NLS 
equation appears as an asymptotic limit for a slowly varying dispersive wave 
envelope propagating in a nonlinear medium. 

In the last years there has been a strong interest in different types of NLS 
equations with cubic nonlinearity whose coefficients depend on the evolution 
variable. As the first example we mention stabilized solitons, which appear in 
multidimensional NLS equations when the nonlinear coefficient is controlled 
appropriately. These objects were predicted to exist in optical applications 
[7], but it was soon realized that a similar phenomenon arises in the context 
of mean field models of Bose-Einstein condensation (BEC) [8,9,10]. As the 
second example, relevant to the present paper, is enhancement of the collapse 
by varying nonlinearity [11]. More references on this field of activity can be 
found in Ref. [12]. 

Another field of application where the coefficients of the NLS model depend on 
the evolution variable is that of dispersion-managed optical solitons. In that 
case it is the second derivative term the one which is modulated along the 
longitudinal direction in an effort to compensate the loss and gain dynamics 
(which is also incorporated into a spatially-dependent nonlinear coefficient) 
(see e.g. the review [13]). The spatial aspects of the dynamics of these disper- 
sion managed solitons are another topic of recent interest [14] where modulated 
coefficients arise. 

Stimulated by Ref. [15], numerous studies of the one-dimensional NLS equa- 
tion with combined temporal modulation of various coefficients, including also 
linear dissipation, have recently been reported [16,17], where the main atten- 
tion was focused on properties of solitary wave solutions. 

The aim of the present paper is to describe several two- and three-dimensional 
models which can be reduced to a standard NLS equation with constant co- 
efficients by means of similarity transformations and describe some of their 
solutions concentrating on localized pulsating pulses, blowing up solutions as 
well as on the solutions which collapse in the case of time independent coeffi- 
cients but are global in the case of varying coefficients. 
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2 Similarity transformations for nonautonomous NLS equations 



2. 1 Model statement and equations for the transformation parameters 



To fix ideas we will denote the evolution variable as t, and concentrate on the 
dynamics of solutions of the NLS equation 



^ = + + gm\ 2 ^ - iim (i) 

For the specific situation of applications of this model in the mean-field theory 
of BEC we have added a potential depending quadratically on the spatial 
variables (r = \x\). We will consider complex solutions of Eq. (1) defined on 
R d , i.e. x E W 1 , satisfying initial data *p(x,0) = tpo(x). 

Let us now make the similarity transformation 



^=wf m Hm T{t) )< (2) 

where t (t), L(t), f(i) and r(t) are scaling functions to be fixed later. 

This transformation is analogous to the so-called lens transformation [18] 
which has been also used in other BEC problems [19,20] and differs (mainly 
because of the extra freedom provided by the definition of the new time r) from 
the transformations which are frequent in the context of dispersion-managed 
solitons [15]. For simplicity and without loss of generality we choose to — 0, 
£(0) = 1, L(0) = 1 and r(0) = 0. 

In this paper we will use the nonlinear rescaling group given by Eq. (2) to 
transform the time dependent problem given by Eq. (1) into simpler, time 
independent ones by choosing appropriately the transformation parameters. 
This transformation has been used in the study of linear Schrodinger equations 
in Ref. [21] but in the opposite way, i. e. to tranform the linear Schrodinger 
equation into a time-dependent problem whose asymptotic behavior provides 
information on the autonomous one via rescaling. 

To determine the scaling functions we require the function u = u(r), r), where 
i] = x/L(t) to solve the NLS equation with the stationary coefficients 



.du 1 l2 . . 

i— = --A v u, +a\u\ u (3) 
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where a = ±1 (the choice of a may imply a redefinition of the sign of g) and 
the initial data for Eq. (3) are given by u(r], 0) = e~ l ^° r ipo{f]), being rj = x and 
fo — 7(0)- Then, inserting Eq. (2) into Eq. (1) leads to the following equations 
for the scaling parameters 



di 

dt 
dj_ 

dt 
dL 

~dt 

dr 

~dt 

a(t) 



a{t)fd£ + ~f(t)e, 


(4a) 


-2a(t)f - in(0, 


(4b) 


2a(t)fL, 


(4c) 


a(t) 
L 2 ' 


(4d) 


git) 
a p ■ 


(4e) 



What we have now in Eq. (3) is the autonomous cubic NLS equation. Thus, 
the transformation given by Eqs. (4) allows one to study the properties of 
solutions of Eq. (1) from those of the simpler, cubic, nonlinear Schrodinger 
equation given by Eq. (3) for which many results are available (see e.g. [18]), 
provided the coefficients are such that they satisfy Eqs. (4). In particular we 
can construct explicit solutions of Eq. (1) from those of Eq. (3). 

An obvious application is to start with stationary solutions of Eq. (3) which 
are known to exist for a = —1 [18] and from them construct breathing and 
blowing up solutions of Eq. (1). In this case, all the time dynamics of the 
solutions, which are self-similar, is contained into Eqs. (4). We will present an 
explicit example of this methodology later. 

Another application, valid specially for the case of two spatial dimensions, is 
to use the method of moments in the framework of Eq. (3) to compute the 
evolution of some relevant integral quantities in the non-autonomous case. 
Finally, explicit soliton and multisoliton solutions of Eq. (3) for d — 1 can be 
used to construct solutions of Eq. (1) with more complex dynamics in the case 
of coefficients varying in time, a fact which has been partially explored in the 
field of dispersion-managed solitons. 

These are only a few illustrations of what can be expected from the application 
of our transformations to specific situations using them as a link between Eqs. 
(1) and Eq. (4). 

In what follows we concentrate on describing the possible dynamics which can 
be obtained from Eqs. (4) with the understanding that the full dynamics in 
the framework of Eq. (1) will be a combination of the dynamics of u given by 
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Eq. (3) with the dynamics of the scaling parameters given by Eq. (4). 



2.2 Solution of the equations for the scaling parameters 



The solution of Eqs. (4a) and (4c) is immediate 



£(t) 



r(t)exp (d f a(t')f(t')dt' 



\ Jo 



(5a) 



Lit) 




(5b) 



where 



T(t) 




(6) 



and thus obviously £(t) = T(t)L(t) d / 2 which, together with Eq. (4e) implies 
that 



This means that the transformations proposed allow us to obtain solutions 
for any of the functions a(t), g(t), and j(t), but once two of them are chosen 
the third one is fixed by (7). Although this change is only valid for a limited 
number of situations we wish to remind that in those remarkable cases the 
information we get relating solutions of the explicitly time-dependent problem 
with to those of the autonomous Eq. (3) is very complete. 

To proceed further we need to solve Eq. (4b) since, using Eqs. (5), it would 
give us the values of L(t),£(t). Unfortunately the general solution of Eq. (4b) 
cannot be written explicitly except for very special cases. 

We will concentrate first on several specific examples which are specially in- 
teresting for applications and will return later to the full variety of behaviors 
which are possible within the framework of Eqs. (4). 



g(t) = aa(t)T(t) 2 L(t) 



d-2 



(7) 
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3 Applications 



3.1 Systems without external potentials Q(t) = 0. 



This situation is typical of problems arising in nonlinear optics (i.e. dispersion- 
managed systems) in which a(t) is a periodic piecewise constant function of 
the propagation variable (which is usually denoted z because of its physical 
interpretation as a propagation distance). 

Let us first consider a dissipationless case when Q(t) = 0, Eq. (4b) can be 
integrated to get 



fit) 2f J^a(t')dt' + l (8) 

Substituting this expression into Eqs. (5) we get the solution in the form of 
quadratures. 

A typical problem in the context of dispersion-managed systems is the ex- 
istence of periodic solutions. Limiting the consideration to the case where 
2/o Jo ot{t')dt' < 1 and using explicit forms of the solutions we get that in 
order to satisfy the periodicity condition L(T) = L(0) the following condition 
must hold 



a{t) dt = 0. (9) 



o l + 2f J*a(t>)dt> 



Using the change of variables (3(t) = fo Jq a(t')dt' this condition can be seen 
to be equivalent to log (1 + 2(3{T)) = and then /3(T) = which leads to 



T 

a(t')dt' = 0. (10) 



o 



In particular, this result implies that under the assumption Q(t) = any 
function a(t) periodic and with zero average will lead to periodic L(t),£(t) 
and f(t). From now on, we assume that T is the smallest period of a(t). 

In a general case using (8) we can integrate Eq. (4d) to obtain 



U 2/oJ*a(f)df + l' 1 ; 
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Subject to the constrain (10) we have that r(T) = and r(t) is bounded, i.e 
< r(t) < t m with T m r M < 0. 



3.2 Systems without dispersion management a(t) = 1 

This situation arises physically in Feschbach-resonance managed Bose-Einstein 
condensates for which g(t) can be continuously controlled but the dispersion 
is fixed a(t) = 1. Obviously, in the case ol = 2 g = aT(t) 2 and sources of 
nontrivial dynamics are the trapping potential and dissipative term depending 
on time. 

In all situations with a — 1, the relevant equation to solve is the Riccati 
equation (4b) which in our case reads 

f t = -v*-\m. (i2) 

For arbitrary dimension d, some general information can be obtained from 
classical tools from the qualitative theory of ordinary differential equations. 
For instance, if f2 is T-periodic and negative, then there exists two (real) 
periodic solutions of Eq. (4b), one of them positive and the other negative. This 
is easily deduced from the classical theory of upper and lower solutions [22]. 
This sign information can be used to get the qualitative asymptotic behaviour 
of £, L from Eq. (5). For instance, if / is chosen positive then 



lim l(t) = lim Lit) = +oo 
whereas if / is negative then 

lim £(t) = lim Lit) = 0. 



(13a) 
(13b) 



Additional information can be provided by writing Q(t) = A + Cl(t), where 
Jq Cl(t)dt = and A is considered as a parameter. Then, a typical Ambrosetti- 
Prodi result holds, that is, there exists a Hopf bifurcation for a critical value 
Ao < 0, such for A < Ao there exist two periodic solutions, while for A > Ao no 
periodic solutions exist. Finally, for the critical value A = Ao, there is exactly 
a periodic solution. 

When Q is a general function (not necessarily T-periodic), the following rea- 
soning is at hand. If x is a real solution of the linear equation 

x + tt(t)x = 0, (14) 
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then it is a simple substitution to check that f(t) = x(t)/2x(t) is a (real) 
solution of Eq. (12). Note that from (5) we get £(t) = [x(t) /x(0)] d/2 , L(t) = 
x(t)/x(0). Eq. (14) is the well-known Hill equation, which has been widely 
studied [23]. If Q is T-periodic, then Eq. (14) is oscillatory (i.e., every solution 
has an infinite number of zeroes going to ±oo) if and only if the mean value 
of Q is positive. In this case, the solution f(t) = x(t)/2x(t) is only defined 
in a finite interval (t ,ti) where t , ti are two consecutive zeroes of x. As a 
consequence, 



lim £(t) = lim L(t) = 0. (15) 
t— >ti t— >ti 



As the first practical output of the above analysis we consider an "explicit" 
example of a pulsating solution in the critical dimension d = 2. Indeed, let us 
consider 



n(t) =m(l -2sn 2 (t,m)), (16) 

where sn(x,m) is a Jacobi elliptic function. A solution of Eq. (14) can be 
constructed explicitly in the form 



x(t) = dn(x, m), (17) 
thus x(t) is a positive periodic function. Let also U be a solution of 



^A V U - U + U 3 = 0, (18) 

with U{rf) ^0 when \q\ — > oo and u(r), r) = U(rf)e lT . Such kind of solutions is 
known to exist [18]. For instance one of them is the so-called "ground state" 
(or Townes soliton) which is a nodeless solution decaying exponentially at 
infinity. Then, the function 



, 1 tt x u\\ ir(t) f . 2 mcn(t,m)sn(t,m)\ 

ib = —— -U — : r, Tit) e lT(t) exp ir l ^— - — — - , (19) 

^ dn(t,m) V dn (^ m ) / V 2dn(t,m) J 

where r(t) = J * dn(t', m)~ 2 dt', is a periodically pulsating solution solution of 
the equation 



= -\ + y (1 - 2 sn 2 (t, m))x^ - (20) 
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Another practical implication, which can be obtained from the above analysis, 
is related to collapse in three-dimensional scenarios for any of the two situa- 
tions described by either Eqs. (13b) or (15) (with the difference, that in the 
first case collapse happen at infinite time while in the second case it occurs at 
finite time). To this end, first we notice that either in a dissipationless case or 
in the case of a periodically varying dissipaion with zero mean value: i.e when 
So ' l(t)dt = 0, if d — 3 and because of Eq. (7) L(t) — > implies g(t) — > 0. 
Moreover we can take as initial data any stationary solution of Eq. (3) (i.e. 
the ground state [18]) when a = — 1. Then we have self-similar collapsing so- 
lutions with decaying nonlinear coefficient given by Eq. (7). In this situation 
what we have is an example of a three-dimensional scenario with a decaying 
nonlinearity for which the initial data collapses as a whole, i.e. an example 
of induction of weak collapse in a three-dimensional scenario. This leads to a 
situation characteristic of two-dimensional scenarios for which the dynamics 
is "less singular" than the usual strong collapse scenarios happening in three- 
dimensional problems with cubic nonlinearities. The smoother behavior of the 
situation is induced by the decay of the nonlinearity near the singularity. 



4 Some general results for sign definite a(t) 

In this subsection we consider we general case under the restriction of definite 
a(t) 7^ for any t. Again, upper and lower solutions together with standard 
arguments from the theory of ordinary differential equations provide some 
useful information. 

If a, n are periodic of the same period and a(t)Q(t) < for all t (that is, 
the respective signs are constant and opposite) then there exists two (real) 
periodic solutions of Eq. (4b), one of them positive and the other negative. 
Also a Hopf bifurcation like in the previous case holds, where the sign of a 
determines the direction of such a bifurcation. 

Other line of procedure is the following. The change in the dependent variable 



V = -7X + «(*)/. 



(21) 



transforms the Riccati equation (4b) into the simpler one 



£ = -V-ift(«>, 



(22) 



with 
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(23) 



This transformation allows to reduce the analysis to the case of constant dis- 
persion a(t) = 1 which was discussed in Sec. 3.2 

Finally, we want to discuss how, given an arbitrary T-periodic coefficient a(t), 
we can design Q(t),g(t) in order to get periodic responses as those obtained in 
the system without external potential. Looking to Eq. (4d), one realizes that 
a necessary condition for the existence of a T-periodic r(t) is that a(t) must 
change its sign. This case cannot be treated with the previous arguments, but 
we can use an inverse procedure. Take L a T-periodic function such that 



r^* = a (24) 

Then r(t) is T-periodic. Once L is fixed, / is determined by Eq. (4c) and 
finally f2 is given by Eq. (4b). Of course, g is always fixed by the restriction 
Eq. (7). 

Finally, if we are not interested in the periodicity of r(t) but only in the rest 
of coefficients, a similar trick is at hand. Given an arbitrary a, take / such 
that a(t)f(t)dt = 0. This guarantees the existence of periodic solutions £, L 
of Eqs. (4a)- (4b) (in fact all the solutions of these equations are periodic) and 
again Q is determined by Eq. (4b). 



4-1 Systems with quasiperiodic coefficients. 



Let us suppose that the coefficients a, Q are periodic with different period 
or, more generally, quasiperiodic functions. Assuming again that a(t) ^ for 
any t, the change given by (21) leads again to the Riccati equation (22) with 
a quasiperiodic coefficient (l(t). Then, the general theory of almost-periodic 
equations [25] can be used in order to prove a Hopf bifurcation as in Subsection 
3.2. The analogous of the mean value for a quasiperiodic function p(y) is the 
generalized mean value 



M lP] = T lim 7^ f P(t)dt- (25) 

Again, we parametrize £l(t) = A + p(t) with M[p] = and then there exists 
A < 0, such that there exist two quasiperiodic solutions for A < Ao and no 
quasiperiodic solutions for A > Ao- In this situation the information on the 
critical value A = A is lost, and the only result remaining is the existence of 
a bounded solution. 
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5 Conclusions 



In this paper we have constructed similarity transformations connecting some 
families of Nonlinear Schrodinger equations with time- varying coefficients with 
the autonomous cubic nonlinear Schrodinger equation. 

These similarity transformations hold when specific conditions linking the 
modulation of the dispersion, potential and nonlinearity are satisfied. In those 
cases they can provide valuable information since they allow to apply all known 
results for the autonomous cubic nonlinear Schrodinger equation to the non- 
autonomous case. 

In particular, using stationary solutions of the autonomous nonlinear Schrodinger 
equation we have discussed how to construct exact breathing solutions to the 
non-autonomous nonlinear Schorodinger equation. 

In addition to discussing the dynamics which can be expected in several cases 
we have provided an specific application in which we explicitly construct 
time dependent coefficients leading to solutions displaying weak collapse in 
three-dimensional scenarios, which is a very unexpected result, i.e. we ob- 
tain a family of collapsing solutions with a behavior which is characteristic 
of two-dimensional scenarios in a three-dimensional setting by taking a time- 
dependent nonlinear coefficient which vanishes at the collapse point. 

Finally let us comment that an extension of the model equation in which only 
p < d of the spatial variables are modulated, i.e. with spatial derivatives of the 
form a(t) J2 P j=i d 2 /dx 2 + J2 < j= p +i d 2 /dx 2 instead of the Laplacian operator, can 
also be treated following the same methodology presented here. However, the 
analysis becomes more complicated and the richness of scenarios is drastically 
reduced in relation to what is found here since the unmodulated derivatives 
impose strict conditions on the possible variations for g(t). 

Our results can find physical applicability in specific examples of mean field 
models of Bose-Einstein condensates and in the field of dispersion-managed 
optical systems. 
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